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The Schwinger-Dyson Equations (SDEs) of matrix models are known to form (half) a Virasoro 
algebra and have become a standard tool to solve matrix models. The algebra generated by SDEs in 
tensor models (for random tensors in a suitable ensemble) is a specific generalization of the Virasoro 
algebra and it is important to show that these new symmetries determine the physical solutions. 
We prove this result for random tensors at large N. Compared to matrix models, tensor models have 
more than a single invariant at each order in the tensor entries and the SDEs make them proliferate. 
However, the specific combinatorics of the dominant observables allows to restrict to linear SDEs 
and we show that they determine a unique physical perturbative solution. This gives a new proof 
. that tensor models are Gaussian at large N, with the covariance being the full 2-point function. 
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^T) ■ INTRODUCTION 

Matrix models have provided a good description of two-dimensional quantum gravity coupled to matter (equiva- 
lently, non-critical strings) [1]. In the scaling limit, they provide access to Liouville gravity coupled to critical (unitary 
' or not) matter with central charge c < 1. While there exist various ways to solve matrix models, the Schwinger-Dyson 
(D . Equations (SDEs, also known as loop equations) have become a standard powerful tool. They enable to probe the 
■ correlators at all orders in the 1/N expansion and the topological expansion [2] has been developed as an intrinsic 
method to solve them. In the double-scaling limit, the SDEs are equivalent to the string equation, which in turn 
corresponds to some integrable hierarchies (depending on the model) [3-7]. 

However the SDEs should not be only considered as a tool, as they probe the key features of matrix models very 
^sO . deeply. In particular, they can be recast in terms of differential operators which generate (half) a Virasoro algebra. 
This naturally leads to the fermion gas formalism which makes the initially hidden, conformal symmetry manifest, 
and consequently clarify the relationship of matrix models to 2d conformal field theories [8-10]. 

It is natural to extend those results to higher dimensions and the proposal was made in the early 90s to use tensors 
instead of matrices (objects with more indices) [ ]. Rank d tensor models indeed generate discrete spacctimcs in 
dimension d and are thus of high interest in different approaches to quantum gravity [14-17], e.g. (causal) dynamical 
triangulations [18-22] and loop quantum gravity [23] (and more generally sit as a natural generalization of matrix 
models wherever the latter are relevant - see [24] for an application to disordered systems). 

A large N limit has been found only quite recently for tensor models [25-27] , but the subsequent developments have 
rapidly expanded. The large N contributions are specific discretizations of the c?-sphere, known as melonic [28], which 
provide an analytic description of the universality class of the Branched Polymer (BP) phase of Euclidean dynamical 
triangulations [30, 31] (previously known from numerical simulations). Moreover, the multi-critical behaviors [29] can 
be interpreted as critical, non-unitary matter [32, 33], just like in 1-matrix models [34]. 

Tensor models not only reproduce the statistical properties of the BP phase, but provide a new way to understand 
dynamical geometries through the SDEs. The algebra they generate was found at large N in [36], and extended at 
all orders in [37]. Geometrically, while the loop equations of matrix models describe the disc amplitude, the tensor 
SDEs describe the ball amplitude. The generators are labeled by boundary triangulations. 

As a first step towards a better understanding of this new symmetry algebra, we prove that the SDEs admit a 
unique physical solution in perturbations. This means that the symmetries completely determine the solution of the 
model, similarly to the conformal symmetry in two dimensions. This unique physical perturbative solution is obviously 
the same as found in [29, 38], by means of scaling arguments and a precise investigation of the combinatorics of the 
leading order contributions in the Feynman expansion. Thus, our method has the advantage that it bypasses the 
Feynman expansion, just like solving the loop equations at large N does not require the knowledge of the planar 
sector in matrix models. 
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The main difficulty compared to matrix models is the proliferation of obscrvables. Indeed, at each order in the 
tensor entries, tensors allow more than a single U (N) invariant (i.e. several boundary triangulations with the same 
number of simplices) and the SDEs are precisely equations on the expectation values of these obscrvables. However 
only a specific subset of observables is relevant at large N and our analysis will consider this family as an input. The 
structure of the relevant observables is such that we can focus on linear equations, in contrast with matrix models. 
This linearity provides an alternative way to [38] to explain why tensor models are Gaussian at large N, the covariancc 
being the full 2-point function 1 . Though the sub-leading corrections to expectation values in the 1/N expansion are 
not known, it is clear that non-linearities eventually come into the game. We will restrict here to the large N limit. 

The organization is as follows. In the Sec. I we briefly review tensor models and the universality theorem. The 
SDEs are derived in the Sec. II, including our fundamental subset of linear equations. Our main result (the unique 
physical perturbative solution) is given in the Sec. Ill and we offer a precise comparison with the loop equations 
of matrix models. The Sec. IV is an attempt to give a global view on the solutions of the SDEs by comparing the 
number of observables to the number of independent equations. This is done in a simple example where the relevant 
observables are shown to be 1-to-l mapped to non-crossing partitions and we conjecture that an infinity of 'initial 
conditions' is required. 



I. A BRIEF REVIEW OF TENSOR MODELS 



Building a tensor model requires first a suitable choice of tensor ensemble, defined by its invariance properties 
(analogously to the Gaussian Unitary Ensemble, or the Gaussian Orthogonal Ensemble in random matrices). A 
natural choice (the only one for which the large N limit is known to exist) is an independent unitary invariance on 
each tensor index. If T ai ... ad are the components of a rank d tensor, at = 1, . . . , N for i = 1, . . . , d, define the following 
transformation, 

bi,—,b d 

where the matrices J/w are independent unitary matrices of size N x N. The complex conjugated tensor T transforms 
with the complex conjugated matrices. We are interested in functions / over such complex tensors that are invariant 
under those unitary transformations, /(T, T) = f(T',T'). They are generated by invariant monomials built in the 
following way [38] : take p copies of T and p copies of T and contract all indices in such a way that an index in the 
i-th position of a T is contracted with an index in the i-ih position on a T. 

Those invariant monomials are conveniently mapped, in a one-to-one fashion, to d-colored bipartite graphs, usually 
referred to as bubbles. Each T is represented by a white vertex and each T by a black vertex, giving, say, p black 
and p white vertices. The indices of each tensor are represented by half-lines labeled by their position, which we call 
color, from 1 to d. When two indices are contracted between a T and a T, they must have the same position hence 
the corresponding half-lines have the same color and one simply joins them together to form a line labeled with that 
color. The invariant monomial obtained from a bubble B can be thought as the 'trace over the bubble'", and we 
denote it B(T, f). 

Bubbles are naturally dual to colored triangulations of d — 1 pseudo-manifolds. The idea is to associate a (d — 1)- 
simplex to each vertex whose half-lines represent the d boundary (d — 2)-simplices of the (d— l)-simplex. Note that the 
(d — 2) simplices inherit the color of the half-lines. Colors further allow to identify all lower-dimensional sub-simplices 
by considering the sub-bubbles with exactly k < d colors. A bubble line between two vertices describe the gluing of 
two simplices along a boundary simplex identified by its color. We refer to [ ] for details. 

Notice that for d = 2, T is a complex matrix and there is only one bubble with 2p vertices: the loop with alternating 
colors 1 and 2, associated to the trace invariant tr(TT^) p . Geometrically, it is dual to a loop with 2p lines. 

Let / be a finite set, {Bi} ieI a set of bubbles and {ti} ie i a set of couplings. A generic action for tensor models is 

S(T,T)=T -f + ^2uBi(T,f), (2) 

iei 

where T • T = ^2 a . T ai ... ad T ai ... ad is the quadratic part (associated with the bubble formed by two vertices connected 
together by d lines). The partition function Z and the free energy F are given by 

exp-N d F = Z = [[dTdT] exp —N^ 1 S(T, f). (3) 



1 Notice that this universality theorem does not prevent critical behaviors as the covariance itself is a non-trivial function of the couplings 
which develops non-analytic singularities [29]. 

2 It was often denoted tr g(T,T) in the literature. 



Figure 1. On the left: an elementary melon with external color 4. On the right: a typical melon on the color 3, with some 
melonic insertions on the internal lines. To get a melonic graph, one can glue melons next to one another and finally close the 
two external lines. 

Such integrals arc usually ' understood as power series in the couplings (perturbed Gaussians). It can be shown that 
F has a 1/N expansion which starts at order 0(1), [29]. 

The natural observables are the bubbles and their expectation values read 

(B(T,f)) = ^ J [dTdf] B(T, T) exp-TV^ 1 S(T,f). (4) 

A distinguished set of bubbles which is of particular importance is the set of melonic bubbles. 

Definition 1. The elementary melon with external color c is defined as the 2-point graph made of two vertices V, V 
connected together by d — 1 lines (having all colors but c), with one external half-line of color c attached to the white 
vertex V and the other, of the same color, attached to the black vertex V (see Fig. 1). A melon is obtained by 
inserting recursively elementary melons on any (internal) line between V and V , starting from the elementary melon 
itself, as in Fig. 1. Melons with the same external color can be joined together so as to get closed connected graphs 
called melonic bubbles. 

Notice that melons in a melonic bubble are precisely the connected, 1-particlc-irrcduciblc, 2-point sub-graphs. A 
melon can also be identified by the two vertices V, V on which the two external lines are attached. This gives a 
canonical way to associate to a black vertex V a white vertex, denoted V. 

The free energy and the Bubble Expectation Values (BEVs) have Feynman expansions onto connected (d + 1)- 
colored graphs (bubbles connected together by propagators which are given a fictitious color) . The way tensor model 
have been solved so far heavily relies on those expansions. The main results have been synthesized in [29]. We recall 
them briefly. 

1. Melonic bubbles. Only melonic bubbles arc relevant at large N. Consequently one can restrict the sum in (2) 
to melonic bubbles only. In the remaining of the paper all bubbles are considered melonic. The 1/iV expansion 
of their BEVs starts at order 0(N), 

(B(T, T)) = N(K B + 0(1/N)). (5) 
Kb is the large N amplitude, given by the theorem below. For a non-melonic bubble, limjv_ i . 0O (B(T, T)) /N = 0. 

2. Gaussian universality. The model is Gaussian at large N, with covariance the full 2-point function. In particular 
the BEVs write 

^(B(T,f))=G^, (6) 

where \B\ is the number of vertices of B and G = (T.T)/N. All dependence on the coupling constants U are 
carried in G. The latter satisfies an algebraic equation (which comes from combining ((>) with a Schwinger-Dyson 
equation) 

1-G-Y^PiUG^ =0, (7) 
where 2pi is the number of vertices of Bi, and with the condition G = 1 when all the couplings go to zero. 



3 About constructive aspects and full summability, see [38]. 

4 A Gaussian theory is defined by the fact that correlators are given as sums over Wick pairings. However in matrix and tensor models, 
not all pairings are equivalent at large N. See the Sec. Ill A for details. 
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cut and remove the vertex 




Figure 2. When the vertex V is cut and removed from the bubble on the left, one obtains a tensorial object, which transforms 
like a T. In this example, d — 3 and (B x V) ai a 2 a 3 = J2b 1 ,b 2 ,b a Ta 1 b 2 b 3 T blb2b3 T bl a 2 a 3 - 

In this paper, we will take the item 1 as granted, as it comes from scaling arguments and amounts to say that we 
have identified the dominant observables. For d = 2 it reduces to a quite trivial statement, as it is equivalent to say 
that the observables are (tT(MM') p ) and that their expectation values start like O(N). In tensor models it however 
becomes a less trivial assertion. We have not found a way to bypass that argument (i.e. derive the dominance of 
the melonic bubbles independently) and we actually find it reasonable to start our study with a given set of relevant 
observables. Then the purpose of the present paper is to find the key equations (6) and (7) by relying only on 
Schwingcr-Dyson equations and without any use of the Feynman expansion of the BEVS onto (d+ l)-colored graphs. 

In addition to the item 1, we will taken as granted the large N factorization 

(B(T, f) B'(T, f)) = (B(T, f )) (B'(T, f )) . (8) 
This is the same assumption that is used in matrix models. It only relies on scaling arguments ". 

II. THE SCHWINGER-DYSON EQUATIONS 

The SDEs and their algebra have been presented in [37]. Since this is not quite standard yet, we re-derive them in 
this section. 



A. Bubble insertions 

The simplest SDE is derived from the identity 

J E [ldTdZ] w ^(T^e-^ 1 <?*+^*<?™)=0. (9) 

^ a 1: ...,a d J 0± ai ... ad 

By taking the derivative explicitly and simplifying by N d , one gets 

1 - (^T ■ f) - Y,PiU (4 Bi(T, f)) = 0. (10) 
iei 

It is a (linear) equation which relates the BEVs together at all orders of the 1/N expansion. It is however not closed 
and one way to close it at large N is to use of the Gaussian universality property (6), which turns (10) into an equation 
on G, namely (7). 

But we have decided not to use the Gaussian universality (and instead to derive it from the SD equations), which 
means we have to write other SD equations. 

Definition 2. Let V be a white vertex in B. The open bubble B \ V, with tensor components (B \ V) ai — a d , is 
obtained by removing the vertex V from the bubble (and the corresponding T ai ... ad in the invariant monomial), so 
that there are d open half-lines carrying the tensor indices a±, . . . , a<i. These half-lines hang out from black vertices, 
hence B \ V transforms like a T. 



It holds gencrically in any theory where the Feynman graphs have faces which bring positive powers of N. The dominant graphs 
contributing to the disconnected part of (B(T, T) B'(T, ?)) have more faces than any connected contribution. 
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An example is given in the Fig. 2. One defines similarly the bubbles B \ V open on a black vertex V, which 
transform like a T. If V S B and V' <E B' , the open bubbles B \ V, B' \ r can be contracted on their free indices 
to get an invariant under (1). We denote it 

(B \V) ■ (B' \ V) = Y, (B^V) ai ... ad (B'^V') ai ... ad . (11) 

ai a d 

This operation has been called bubble gluing in [37]. 

Open bubbles appear naturally in the derivatives of a bubble invariant, 

§^ = £ (B x V) ai ... ad , d -§^± =X>^ F) 0l ... 0d . (12) 

^oj-a,, yeB 01 ai ... ad - gB 

What happens when an open bubble is opened a second time? That gives an object with more indices (and possibly 
disconnected as a graph, i.e. which factorizes as the product of two open bubbles). This can be repeated several 
times. Conversely, one can take a bubble opened several times and sum over some indices as long as one contracts an 
index of a T with one of a T (to satisfy (1)). This operation has been called bubble contraction in [ V]. 

Open bubbles are covariant objects that can be used as insertions to generalize (9). SDEs are thus labeled by open 
bubbles and come from the identity 



ai,...,a d J oi—a,j 



(13) 



We have to distinguish three types of contributions. 



1. The simplest one is when the derivative acts on the quadratic part of the action. Indeed d(T ■ T)/dT ai ... ad = 
f ai ... ad and by definition, J2 ai ,....a d ( B x V) ai -a d T ai ... ai = B(T, T). Therefore 

Z a u ...,a d J CJ 0l ...„ d 

2. Another contribution comes from taking the derivative of the bubble terms of the action. Using (12), the 
derivative of a bubble Bi acts on its pi white vertices Vi and produces for each of them the open bubble Bi \ Vi 
which transforms like a T. Hence 

E fidTdf](B^V) ai ... ad YU^^e-^- 1 ^^ 3 ^ 

ai,...,a d i£l oi—Od 

i€l VitBi \a 1 ,...,a d / 



= -N d - 1 Y* l E ((B^V)-(Bi^Vi) 

iel V z EBi 

Here the sum over Vi runs over the white vertices of Bi. When B and Bi are melonic, all the resulting bubbles 
(B \ V) ■ (Bi \ Vi) are melonic, with \B\ + 2pi — 2 vertices, and the scaling is the same as in (14). 

3. The last contribution is the derivative of the open bubble B \ V. It produces a sum over its \B\/2 white 
vertices, where for each term the white vertex is removed. That opens the bubble a second time, and the 
resulting open half-lines carry the indices ax, ■ ■ ■ , &d- They transform like a T, and performing the sum over the 
indices connects these half-lines with those of B \ V, producing an invariant which we denote B \ V \ V (this 
is a bubble contraction). We can write 

\ Y [[dTdf] ^_lYh^i e -N d - 1 (T-T +I : ieI uB i (T,f)) = ^2 (B\Ky'(r,T)). (16) 

Z a 1 ,...,a d J °T ai ... ad yieB 

Each B \ V \ V is a graph with two vertices less than B, and is typically a disconnected set of bubbles (see 
below) . 
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Figure 3. On the left: the bubble B with the canonical pair V, V and 2-point insertions Mi, . . . , Md, for d — 3. On the right: 
the open bubble where the striped circles are melonic insertions (which make the structure of Mi , . . . , Md more explicit) . 

Cutting and removing one white vertex leads to disconnected components. To reach the maximal scaling N d , d disconnected 
pieces are required. The only way to get them is to cut and remove the vertex V canonically associated to V. 



B. The fundamental large N equations 

One can write an exact equation, which holds at all orders in the 1/N expansion, by summing (1G), (14) and (15) 
together. But there is a simplification at large N, as we have to keep only the melonic contributions in (16) which 
scale like (14) and (15), i.e. 0(N d ). It turns out there is only one such contribution, and this can be proved as follows. 

For a melonic B, the vertex V is part of a canonical pair (V, V) and one can always draw B like on the left of 
the Fig. 3. The lines of color 1, . . . , d carry 2-point insertions Mi, . . . , Md that we first assume to be non-trivial. 
When V is removed, one gets the open bubble on the right of the Fig. 3, where Mi, . . . , Md have been expanded to 
describe all typical white vertices. To get B \ V \ V, a white vertex V has to be removed. Each disconnected piece 
in B \ V \ V scales like N, so that to reach the scaling 0(N d ), B \ V \ V must contain d disconnected pieces. 
It can be checked explicitly using the figure that there is only one way to get d disconnected pieces, which happens 
for V' — V, the vertex canonically associated with V. The bubbles which result from the contraction are simply the 
2-point sub-graphs Mi , . . . , Md which are closed by joining their two external lines (abusing the notation, we still 
denote the corresponding bubbles M c , for c = 1, . . . , d). Using the large N factorization (8), we are finally led to 

d 

\[(M C (T,T)) - N*- 1 (B{T,T)) - N^Y^U ((-B \ V) ■ (Bi \ Vi)) = 0. (17) 

c=l ieJ Vi£Bi 

If some of the M c s are trivial, i.e. they have no vertices and just consist of lines of color c, closing them produces 
loops with no vertices and free sums on the indices a c , each of which simply producing a factor N. 
The main difficulties with this set of equations are: 

• the non-linearities, which come from taking the derivative of the open bubble, 

• the proliferation of observables: from the bubble B, the equation generates all the possible gluings of B \ V 
with the bubbles Bi contained in the action. This feature contrasts with the loop equations of matrix models 
by bringing an additional combinatorial difficulty. 

The melonic family of bubbles is constructed by recursive insertions of the elementary melon on any line. It turns 
out that this is the reason why it is going to be sufficient to focus on the SDEs for B \ V where V is the black vertex 
of an elementary melon. The equations are given in the following lemma. 

Lemma 1. Linear equations for elementary melons. A melonic bubble B always has at least one elementary melon 
M, with external vertices V, V, and without loss of generality, external color d. The SDE for B \ V involves the 
2-point subgraphs M c , c = 1, . . . , d — 1, which are just closed lines with no vertices, hence contributing as N 11 ^ 1 . 
Define the bubble B \ Md as B with the elementary melon Md replaced by a line of color d. Then (17) becomes 

(B n Mi) - (B) -J^ti Yl ^ B x ^) ' ( Bi x V ^>) = °' ( 18 ) 

These are linear equations which describe the creation/annihilation of an elementary melon Md- The fact that this 
is a complete set at leading order and that they are linear will quite directly lead to the universal property (6). 

An interesting, direct application of the above Lemma is given below (but it will not be used in the proof of the 
universality result). 
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Application 1. Equalities between sums of BEVs. Let G be a 2-point (melonic) graph, say, with external color 1. 
Then, 



E** E ( 2( ^1^ V J 1 §^)) = E^ E C^T^C^ 2 )- (19) 

The proof goes as follows. Consider the bubble B obtained by closing the graph G with two additional elementary 
melons next to each other on the color 1. We call V and V the two white vertices of these melons, and compare 
the SDEs for B \ V and B \ V'. Both contain the expectation value of B. Since in the two cases B is opened on 
elementary melons, the lemma 1 tells us that the SDEs generate a single bubble with two less vertices. This is the 
same bubble in both cases: B with one elementary melon replaced by a line of color 1. Therefore, the only difference 
comes from the gluings of B with the bubbles Bi of the action. The proof ends by substracting one equation to the 
other. 

This application becomes more interesting in special cases. For instance if the bubbles in the action have the same 
number of vertices, then the equation compares the BEVs of bubbles which all have the same number of vertices. 
It is even better for instance if there are some symmetries on the bubbles which reduce the number of terms of the 
equation. 



III. THE GENERIC 1-TENSOR MODEL 
A. The Gaussian model 

Before analyzing the generic model, it is useful to understand the behavior of the Gaussian one. For a covariance 
<7, the action is 

S(T,f) = -T-f. (20) 

The SDEs of the Lemma 1 are particularly simple, 

(B) =g (Bn M) , (21) 

for any elementary melon M in B. As any melonic bubble comes from inserting an elementary melon on a smaller 
bubble, we get with the initial condition (1) = 1 that 

(B)=Ng», (22) 

where p is the half number of vertices. This result has a combinatorial interpretation. Indeed a Gaussian measure is 
standardly defined by the fact that the expectation values are sums over Wick pairings weighted by the covariance. 
However, in matrix and tensor models not all Wick contractions have the same scaling with N, which means that 
some are suppressed at large N. As well-known, in the Gaussian matrix model, only planar contractions survive at 
large N (and are counted by the Catalan numbers). For the Gaussian tensor model, the BEV (22) shows that there 
is a single Wick pairing which survives the large N limit. It has been shown to be the pairing which connects the 
canonical pairs (T, T) of the melons in B 6 . 

B. Solving the generic 1-tensor model 

In this section / is a finite set, {Bi}i e j a set of bubbles with associated couplings {ti},^/. We denote the half-number 
of vertices of each bubble by pi = \Bi\/2. We will need the following Proposition. 



In other words, there is a single way to build a melonic graph with (d + 1) colors from a melonic graph with d colors and this is the only 
way to maximize the number of faces. This is the combinatorial reasoning used in [38]. 
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Proposition 1. Let (F n ({ti})) n ^n be a sequence of functions such that each F n has a power series expansion in each 
ti. Let E be a map which associates to each bubble B a series Esiiti}) in the couplings {ti} such that the empty 
bubble is mapped to the constant 1. Assume that for any vertex V of an elementary melon M in a bubble B with 2p 
vertices, E satisfies the equation 

Eb-^m - E B 2J e (b^v)-(b,^v') = f p- ( 23 ) 

Then the evaluation E on two bubbles with the same number of vertices gives the same function, i.e. E only depends 
on the number of vertices of the bubbles B and not their specific structure. 

Proof. We proceed by induction on the number of elements in /. 

• When I is empty, Eb^m — Eb only depends on the number of vertices of B, hence the result follows from 
a trivial recursion on the number of vertices of the bubbles (as any bubble can be obtained by inserting an 
elementary melon on a smaller bubble). 

• Assume the lemma holds for a particular set / and consider an additional bubble Bo, with 2po vertices and 
coupling t . We expand Eb as 

E B (t ,{t i }) = ^EP({t i })- (24) 

n>0 

The function F p has a similar expansion as F p = Yln^o^v • The equation satisfied by E becomes at order 
n > 1 

p (n) _ p (n) _ , F (n) _ p(n) , / 9 r\ 

ie/ v'eBi v&Bo 

and 

p(o) p (o) v^f - ^fi^ 

B-^M B 2^ 1 £ '(B\V)*(B,\V') ~ P ' ^ ' 

ie/ y'GB; 

At this stage, we would like to apply the lemma for the set / as our induction hypothesis. However the right 
hand side of (25) is a priori not a sequence since it seems that it depends on gluings of bubbles. However, we 
now show by induction on the order of the perturbation that it only depends on the number of vertices of B 
and that the lemma for the set / can be applied to (25). 

— At n = 0, the lemma can be clearly applied on (2(i) with the set /. As a result (Eg\{ti})s is actually a 
sequence of functions since its sole dependence on B is through the number of vertices. 

— Assume this holds at order n — 1, for n > 1, we apply the lemma to (25) for the set I and find that the 

fn) 

dependence of E B on B is just on the number of vertices of B. 

Therefore, the full expansion (24) of Eb only probes bubbles through their number of vertices, which proves the 
desired property for the set supplemented with Bq. 

□ 

We now apply the above proposition to the tensor model with the generic action 

S(T,f) = T-f + ^2u tr Bi (T,r). 

iei 

One sets E B = {B(T, f )), F = and the large N BEVs satify the equation (23) (the Lemma 1). 

Corollary 1. The bubble dependence of the BEVs is just through their number of vertices. For any bubble with 2p 
vertices, denote (B(T,T)) = N G p . Then the linear SD equations reduce to 

G n - G n+1 ~Y,Pi *i G n+ P , = 0- (27) 
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Let us re-organize the set of bubbles {Bi}i £ i according to the number of vertices of the bubbles. Set I = U pe pJ p , 
where P is a hnitc set of integers greater or equal to 2, such that {Bi} ie j p is the subset which contains the bubbles 
with 2p vertices. We define the coupling at 2p vertices as t p = J2iei p Then the recursion reads 

Gn ~ Gn+l — *^^ptpG n + p = 0. (28) 
pEP 

This linear recursion can be solved in two (equivalent) ways. 

1. Solving the recursion. The characteristic polynomial is 

Pt P X P + X - 1 = p*t p , (X - G) H(X - G (a) ). (29) 

p£P a 

Here p* = SupP is the maximal number of vertices among the bubbles in the action. When the couplings {t p } 
go to zero, there is a single root X = 1. We have denoted G the root which goes to 1 when the couplings go 
to and we have isolated it on purpose as it is the physical root. Assuming that the P roots G, (Gr a \) of this 
polynomial arc distinct, it comes for any B with 2p vertices 

(±B(T,f))=cGr + J2c«GU> ( 30 ) 

for some constants c and c a . The physical interpretation is clear from the analysis of the Gaussian model in the 
Sec. Ill A: each geometric contribution G^ a \, G p , is a large N, Gaussian channel with covariance G( Q ),G. But 
only one of them goes to 1 when the couplings go to and it is the physical Gaussian channel with covariance 
G. Therefore (B) = N G n , and G is the full 2-point function. Note that this argument avoids the discussion of 
the initial conditions which are in principle necessary to solve the recursion. 

2. With the resolvent. Natural objects of matrix models are the resolvent and its associated eigenvalue density. 
While tensors have neither natural multiplication like matrices, nor eigenvalues, the fact that the large N BEVs 
only depend on an integer (the number of vertices) enables to define a natural analog to the resolvent 

u{z) = Y,^ P ~ l G p , (31) 

p>0 

In matrix models it turns the SDEs into an algebraic equation on lo(z). We can proceed similarly here, to get 

G Q + YJ P p=iP t pGp-i +YJk=i ^YTp^k+iPtpGp-i-k 
= ^ : — — 7 , (32 

where G& for fc = 0, . . . ,p» — 1 are the required initial conditions. To avoid setting the initial conditions by 
hand, we need to identify the physical Gaussian of covariance G (the one which goes to 1 when the couplings 
go to zero) in the language of the resolvent. This is done by analyzing the structure of its singularities. 

We observe that the denominator is the characteristic polynomial (29), whose roots are G, (G( Q )), hence 

G + YZ*=2P t pG P -i +YX=i ^YZ'^k+iPtpGp-i-k 

w O) = : — / ^ n 7 — s ■ (33) 

jj* t p , (z-G)[[ a {z-G {a) ) 

Since the numerator is a polynomial of degree p» — 1, the resolvent has at most p* poles, z = G, z = G( Q ), 
and at least one pole. As only G is physically relevant, we make a 1-pole hypothesis: the polynomial in the 
numerator must remove the p* — 1 non-physical singularities located on (G( a )). This fixes the initial conditions 
Gk, k = 0, . . . ,p» — 1, up to a global scale. The latter is fixed by the trivial condition Go = (1) = 1. This 
reasoning leads to the final form of the resolvent 

w(z) = j^. (34) 

The associated 'eigenvalue distribution' p(X) is obtained as usual by taking the discontinuity of the resolvent 
across the real line, 

p(\) = 5(\-G). (35) 
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Figure 4. Such bubbles can be re-written as loops observables in a matrix form tr(TT t ) p , with p = 3 in the picture, for a 
matrix T a A where a = 1, . . . , TV is the index on the color 1, and A ranging from 1 to TV d_1 describes all the colors 2, . . . , d, with 
d = 4 here. 

C. Comparison with matrix models 

The above resolvent and 'eigenvalue distribution' describe a system of non-interacting particles, all falling in the 
same potential well at z = G. This is in contrast with matrix models, but it is actually what would happen in a matrix 
model if the famous Vandermonde contribution, which acts like a Coulomb gas repulsion between the eigenvalues, 
could be removed. Such a matrix model has actually been constructed in [35]. Further, the universality property (6) 
ensures that all physical quantities of the generic 1-tensor model can be evaluated with this model. 

It is a matrix model for 'very rectangular' matrices, of size ./V x TV that is built as a tensor model with very 
specific bubbles. The bubbles entering the action arc chosen such that the indices a2, . . . , a<j of T aiCl2 ... ad are always 
contracted all at the same time with the corresponding indices of a single T , as in the Fig. 4. The associated bubbles 
look like loops with some lines of color 1 and the others being 'fat lines' with colors 2, . . . , d. These invariants can be 
written in a matrix form. Introduce a 'fat index' of size TV d_1 , denoted by a capital letter, A = (a 2 , . . . , a<j) and write 
the tensor like a matrix (T a) of size TV x TV^ -1 and its complex conjugate T a A = (T^)ao,- The corresponding action 
is 

p, 

S(T, Tt) = tr(TTt) + ^ t p ti(TT^. (36) 

p=2 

This model was analyzed in details in [35]. Let us reproduce briefly the points that are most relevant to us. Upon 
introducing the eigenvalues (Xi)i—i t ,,, t N of TT\ the partition function reads 

r N ( N n / N 2 \\ 

Z = JlldXi cxp-TV^ 1 E^ + ^IE^'- AfI>l Ai - A 'l) I' (37) 

with 



N d-i N 

^i— 1 i<j 



V(X) = A + J]tpA p -In A. (38) 

p=2 

The two In A terms are due to the fact we work with a rectangular matrix (they cancel each other for a square 
matrix). One recognizes X^<j m — Xj | as the famous Vandermonde contribution. It is of order N so when d = 2 
it balances the potential J^. V(Xi). It is the usual matrix model philosophy: the Vandermonde term acts like a 2d 
Coulomb repelling which prevents the eigenvalues to all fall in the minimum of the potential and hence spreads the 
eigenvalue distribution. However in our case, the Vandermonde is rescaled by N/N d ^ 1 and is therefore suppressed as 
soon as d > 2. Consequently the saddle point equation decouples the eigenvalues and reads 

V'(X) = j(jr P t p Xi' + X-l) =0. (39) 

p=2 

One recognizes the characteristic polynomial (29) of the BEV recursion. Thus, the large TV Gaussian channels with 
covariance G, (G( a )) identified in the tensor model analysis naturally appear here as the solutions of the saddle point 
approximation. 

Let us now show how the loop equations describe the above phenomenon. We start with the identity 

/ [dTdT*] -7—— — f \(TT*) n T\ exp-N d - 1 S(T,T*) \ = 0. (40) 

J OT a A \ l ■ ' i aA J 
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Next we evaluate explicitly the derivatives. Using the large N factorization (tr(TTt) fc tr(TTt)?) = (ti-(TT^) k )(tr{TT^)P), 
one gets 

P n—l 

(Itr(TTt)«)-(itr(TTt)^)-^^ p (Itr(TTt)^} + ] J- I £(I tr(TT^) (I tr(2Tt)-*) = 0. (41) 

All quantities into brackets ( } are of order 0(1), hence for d > 2 the non-linear terms are suppressed. In other words, 
the presence of the Vandermondc contribution to the saddle point equation translates into non-linearities in the SD 
equations. For d > 2 the natural scaling of tensor models removes those non-linearities and the linear equations of 
the Lemma 1 are recovered. 



IV. ABOUT OBSERVABLES AND OTHER SOLUTIONS OF THE SCHWINGER-DYSON EQUATIONS 

We proved in the previous section that there is a unique physical perturbative solution to the SDEs at large N. 
However, nothing has been said about the full set of solutions to the equations. We propose some preliminary analysis 
in this Section, based on a specific, simple example, and conjecture that the set of solutions is determined by an 
infinity of 'initial' conditions. 

We consider two bubbles Bi,B 2 , both with four vertices but different color labels, and two different coupling 
constants t\,t2, 



S(T,T) = T-T + t! 



3 +t 2 



(42) 



Bi(T,T) 



B 2 (T,T) 



Note that the bubbles above have been drawn for d = 3. To get them for generic d > 3, it is sufficient to add lines 
with colors 4, . . . , d whenever there is already a line of color 3, i.e. between all canonical pairs (V, V). The model gives 
a special role to the colors 1 and 2. As a consequence, it is natural to define a specific sub-set of melonic bubbles. 

Definition 3. A bubble is said to have melons on the colors 1 and 2 only if it can be built from the 2-vertex bubble 
by recursive insertions of elementary melons on the colors 1 and 2 only. Let B\ 2 (p) be the set of such bubbles with 
exactly 2p vertices, \Bi 2 (p)\ the corresponding number of bubbles, and Bi 2 = U pe fj£>i2(p). 

To write the SDEs, we need to know the open bubbles obtained from Bi and B 2 . The two black vertices of Bi 
(and B 2 ) give the same open bubble B\ \ V% (and Bi \ V\), depicted in Fig. 2. Therefore, for any choice of B, the 
SDE generates two larger graphs, one with an additional elementary melon on a line of color 1 on B and the other 
with an additional elementary melon on a line of color 2. This implies that if one writes a SDE for an open bubble 
B s V which has melons of external colors 1 and 2 only, the equation only generates other bubbles in £> 12 . We can 
write the SDEs graphically, 




0. (43) 



Quite clearly, these equations actually generate the whole set Bi 2 . 

We will first count the number of obscrvablcs generated by B\ and B 2 from the trivial bubble, and then compare 
with the number of independent SDEs, for graphs with up to fourteen vertices. 



A. Melonic bubbles and non-crossing partitions 

As mentioned above, the sets B\ 2 {p) are actually independent of d > 3. If d > 3, then for any line of color 3 between 
two canonically associated vertices V, V, there are also lines of colors 4, . . . , d between V and V, and this situation 
accounts for all lines of colors 3, . . . , d. By removing the lines of colors 4, . . . , d, one obtains an element of i? 12 for 
d = 3. Reciprocally, given an element of B\ 2 for d — 3, one can add lines of colors 4, . . . , d between all pairs (V, V) to 
go back to the generic case. Therefore we focus on d = 3 in the following. 
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4 4b 4a 

Figure 5. One the left: a NCP on 6 elements with labeled vertices. On the right: an element of B^ip = 6). A melonic bubble 
is obtained by putting back the color 2 on each arc between n a and rib, and the color 1 on the arcs between rib and (n + l) a , 
and removing the vertex labels. 

One can represent an element of 612 (p) as a bipartite, planar contraction among 2p elements. Indeed, one first 
places the 2p vertices of the bubble on a circle, alternating black and white vertices. Using the clockwise convention, 
we put the color 1 on the p arcs which go from a black vertex to a white vertex, and the color 2 on the p other arcs. 
Each black vertex must then be connected to a white vertex by a line of color 3, the pairings being allowed only if 
they result in a planar graph. 

The vertices of the bubbles in B12 are not labeled. We introduce B^ as the set of bubbles whose melons have 
external colors 1 and 2 only, but now with labeled vertices. We label them (l a , lh, . . . ,p a ,Pb) going clockwise around 
the circle, where the vertices (n a ) are the white ones and (rib) the black ones, as shown on the right of the Fig. 5. 
Then B12 can be described as the equivalence classes of these labeled planar contractions under the action of the 
rotations (n a t— s> (n + 1) , rib H» (n + l)b) (with periodic boundary conditions), and reflections. 

Let us now introduce the Non-Crossing Partitions (NCPs) of [p] — {1, . . . ,p} in a graphical way. A partition it 
can be pictured by putting the elements of [p] on a circle, going from 1 to p clockwise, and adding links between 
the elements of each sub-set as follows. If n, m form a sub-set {n, m} with two elements, a line is drawn to connect 
them. When a sub-set has at least three elements {k, I, m, . . . }, we draw the unique convex polygon whose vertices 

are k,l,m, Finally singlets {n} are identified as isolated vertices on the circle. The set NCP(p) of non-crossing 

partitions is the set of partitions 7r of [p] with no crossing inside the circle. An example is displayed on the left of the 
Fig. 5. 

We will describe a bijection between B^ip) and NCP(p) and then mod out rotations and reflections. First we 
present the map from NCP(p) to B^ip)- We start by splitting each vertex n = 1, . . . ,p into two, denoted n a ,rib, 
clockwise ordered, and label the arcs between rib and (n + l) a with the color 1 and the arcs between n a and rib with 
the color 2. We can distinguish three types of elements in ir £ NCP(p): 

• for each singlet {n} € 7r, we draw a line (with color 3) between n a and rib (this creates an elementary melon on 
the color 1), 

• for {n, m} G tt with exactly two elements, we draw a line between n a and m^, and a second line between rib and 
m a (notice that they do not cross), 

• for a polygon corresponding to {m, «2, «3, ■ ■ ■ } £ tt, such that the links are between and rik+i, we draw lines 
from (rik)b to (n k +i) a - 

To get the inverse map, one simply reverses the above three rules, which account for all possible patterns inside the 
bubbles. The correspondance is detailed graphically in the Fig. G. Then, it is easy to see using this map that the 
rotations (n a i-> (n + l) a ,?i,b M> (n + l)b) are mapped to usual rotations n rt n + 1 on NCP(p), and reflections to 
reflections. As a consequence, the number of bubbles |i3i2 (p) | is the number of non-crossing partitions of [p] up to 
rotations and reflections (i.e. dihedral classes). They have been studied in [39] and are known as Alll275(p). The 
first values, from p = 1 up to p = 7, are 1, 2, 3, 6, 10, 24, 49. 



B. Independent SD equations 

Now that we have a good control on the family $12, we come back to the SDEs. The key question is: do the SDEs 
determine the BEVs up to a finite number of initial conditions? While we have not been able to give a rigorous 
answer, we conjecture that the answer is no. It is based on the explicit analysis of the equations on £>i2(p) up to 
p = 7. 
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Each SDE relates the BEVs of two bubbles with 2(p + 2) vertices to the BEVs of bubbles which have less vertices. 
From the point of view of the two bubbles in Bi2(p + 2), the equations form a linear system. We will be interested in 
its rank, to compare it in particular with the number of bubbles \B\i(j> + 2)|. Our conjecture is that the rank of the 
system is precisely \B\i(j> + 2)| — 1, for all p > 0. To support it, we first write down the equations for a few values of 
p. For p ~ 0, 



2ti 3 2 2 3 + 2t 2 3 1 1 3 



(44) 



There are two bubbles in £>i2(2), but only one equation. For p = 1, there are two independent equations on three 
bubbles, 




(45) 



For p = 2, there is again one equation less than the number of bubbles |£>i2(4)| = 6 (the notation 1 <-> 2 means that 
another equation is obtained by exchanging the colors 1 and 2, and t\ with £2)1 




and 1 <-> 2, 



and 1 O 2, (46) 



However, when p = 3 is reached, we get 12 equations, which is more than the number of bubbles on ten vertices 
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|Bi 2 (5)| = 10, 




and 1h2, 



and 1 O 2, 



and 1h2, 







<] 


1 











and 1 -o- 2, 



and 1h2, 



and 1h2. 
(47) 

But the system is of rank 9 = |Bi2(5)| — 1 only. A consequence is that it is possible to extract 3 equations which 
do not involve the bubbles of Sia(5) at all and appear as constraints on the BEVs of smaller bubbles. Then they 
might supplement the equations (44), (45) and (46) and determine some of the previously left undetermined BEVs. 
However, this is not the case. We have checked that these constraints are trivially satisfied if (44), (45) and (46) hold. 

We will here refrain ourselves from writing the 35 equations for p = 4 and the 102 equations for p = 5 (though quite 
tedious, it is actually really straightforward to write them in this model since B\ : B2 only add elementary melons). 
In those two cases, the rank is again \B\2{p + 2)| — 1. 

This inspection thus suggests that one BEV must be specified at each order in the number of vertices to determine 
all the others, and that the set of solutions to the SDEs is parametrized by an infinity of 'initial conditions'. To see 
how that could work in practice, we propose to re-organize (some) SDEs according the number of melons on the color 
2. 

When t2 = 0, the model can reformulated as a matrix model like in the Section III C with a simple 1-bubblc 
potential. Then the loop observables with melons on the color 1 obey a special case of the recursion (28). Therefore, 
it can be instructive to understand how a non-zero coupling ti affects this recursion. Denote G n the expectation value 
of the 'loop' bubble with exactly n elementary melons on the color 1. At t2 = 0, for any n > 



G„| t2=0 ~ G„+i| t2=0 — 2ti G„ + 2|4 2= o — 0. 



(48) 



We denote G n ^\ the expectation value of the bubble with n melons on the color 1 plus one elementary melon on the 
color 2 (hence n — 1 elementary melons on the color 1). This is the term generated for £2 7^ 0, 



G n — G n +i — 2t\ G n _|_2 — 2£2 




= 2<2 Gji+i,!. 



(49) 



Let us now focus on the family [G n 1). The bubble corresponding to G n 1 has precisely n — 1 elementary melons on 
the color 1, and another melon on the color 1 whose internal line of color 2 carries an elementary melon. By opening 
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the graph on an elementary melon of external color 1, we see that at ti = it satisfies the recursion (48) for n > 1, 
the reason being that the graphs simply differ by the number of elementary melons on the color 1 , 

G n ,i\t 2 =o ~ G„ + i4| t2=0 — 2ti G n+ 2,i| t2= o = 0. (50) 

When f 2 7^ 0; the recursion picks up an additional term, involving graphs which have two elementary melons of 
external color 2. The BEVs of these graphs (for arbitrary £2) only depend on the number of vertices. Indeed, this is 
a consequence of the Application 1 in the present model. The contribution of the bubble B\ being the same on both 
sides of (19), it gives the following equality 7 

n-k 




We denote the corresponding BEV G n .2, for n > 2 (they have exactly n melons on the color 1 and two on the color 
2), so that 

G n ,i — — 2ti G„+2,i = 2*2 Gn+1,2- (52) 

These results are similarly extended to the family of graphs (G n ,m)n>m obtained by choosing m elementary melons 
on the graph G n and adding an elementary melon on the color 2 on each of them. The equality (51) extends to this 
family, stating that their BEVs are independent of the location of the m elementary melons on the color 2 (as long as 
each of them is added on a different elementary melon with external color 1). This family satisfies a double recursion, 

— 2ti G n +2,m = 2*2 Gn+l,m+l- (53) 

This double recursion does not take into account all SDEs on Bi2- However, it is consistent with our conjecture. 
It is clear it needs an infinite number of initial conditions, namely the family (G n ), which can be freely chosen. It 
then determines the family (G„.i) via (49), which determines in turn the family G n ,2 via (52), and so on. Again, the 
conjecture is then that all other SDEs either determine some BEVs outside the family (G„ )m ), or are redundant. 



V. CONCLUSION 



Our main result is that the large N SDEs of tensor models admit a unique physical, perturbativc solution. This 
provides an alternative proof of the universality theorem (6) with less combinatorics than in the original proof [38]. If 
we think of the SDEs as differential operators acting on the partition function and generating a specific generalization 
of the Virasoro algebra [36, 37], our result actually means that these new symmetries completely determine the 
physical solution. While that was a necessary step, the full representation theory of this new algebra is still to be 
understood. Interestingly, the representation carried by the partition function of tensor models is such that it reduces 
the algebra to a Virasoro algebra, as already noticed in [29] . 

Beyond the large TV limit, we expect the SDEs to give access to sub-leading orders, and possibly to a double-scaling 
limit (where N goes to infinity and the couplings go to their critical values simultaneously). We already know that it 
is possible for the tensor model used in the Sec. Ill C since it is actually a 'very rectangular' matrix model which has 
been shown to possess a double-scaling limit [35] (and the resolvent goes like lo{z) = l/(z — G) + a/N d ~ 2 (z — G) 2 + 
b/N d ~ 2 (z — G) 3 + • • • ). This project is quite interesting since at sub-leading orders the SDEs start to distinguish 
bubbles that have the same number of vertices. It would be useful to find a suitable generating function for bubbles, 
such as the one proposed in [40] which takes colors into account (but it is not clear whether that one is directly 
relevant). Another difficult part is to identify the non-melonic bubbles whose first non-zero contribution starts at a 
given sub-leading order. 

Using those ideas, we hope that will be possible in the future to recast tensor models in a way which puts the 
emphasis on the new symmetries encoded by the SDEs, similarly to the fermion gas formalism of matrix models 
which makes the conformal field theory content clearer. 



7 This proof only works for t2 ^ 0. When t'2 = 0, an independent proof is easily obtained by showing that all these graphs satisfy the 
same recursion, which is actually nothing but (48). 
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